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Abstract 

Gas-flow  dynamics  in  internal  gas  manifolds  of  a  fuel-cell  stack  are  analyzed  to  investigate  overall  pressure  variation  and  flow  distribution. 
Different  gas-flow  patterns  are  considered  in  this  analysis.  Gas-flow  through  gas  channels  of  each  cell  is  modeled  by  means  of  Darcy’s  law 
where  permeability  should  be  determined  on  an  experimental  basis.  Gas-flow  in  manifolds  is  modeled  from  the  macroscopic  mechanical 
energy  balance  with  pressure-loss  by  wall  friction  and  geometrical  effects.  A  systematic  algorithm  to  solve  the  proposed  flow  model  is 
suggested  to  calculate  pressure  and  flow  distribution  in  fuel-cell  stacks.  Calculation  is  done  for  a  100-cell  molten  carbonate  fuel-cell  stack  with 
internal  manifolds.  The  results  show  that  the  pressure-loss  by  wall  friction  is  negligible  compared  with  the  pressure  recovery  in  inlet  manifolds 
or  loss  in  outlet  manifolds  due  to  mass  dividing  or  combining  flow  at  manifold-cell  junctions.  A  more  significant  effect  on  manifold  pressure 
possibly  arises  from  the  geometrical  manifold  structure  which  depends  on  the  manifold  size  and  shape.  The  geometrical  effect  is  approximated 
from  pressure-loss  coefficients  of  several  types  of  fittings  and  valves.  The  overall  pressure  and  flow  distribution  is  significantly  affected  by  the 
value  of  the  geometrical  pressure-loss  coefficient.  It  is  also  found  that  the  flow  in  manifolds  is  mostly  turbulent  in  the  100-cell  stack  and  this 
way  result  in  an  uneven  flow  distribution  when  the  stack  manifold  is  incorrectly,  designed. 

©  2003  Elsevier  Science  B.V.  All  rights  reserved. 
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1.  Introduction 

One  of  the  engineering  problems  encountered  in  the 
design  of  fuel-cell  stacks  is  the  configuration  of  gas-flow 
manifolds  and  the  overall  pattern  of  gas-flow  that  would  give 
uniform  flow  distribution  and  stable  cell  operation.  The 
geometrical  structure  and  size  of  the  gas  manifolds  are 
the  important  parameters  to  be  considered. 

For  planar-type  fuel-cell  stacks,  it  is  first  necessary  to 
decide  whether  internal  or  external  manifolds  are  to  be  used. 
External  manifolds  are  simpler  and  less  costly  than  internal 
manifolds,  but  there  is  a  major  drawback  of  gas  leakage  and 
sealing  [1],  External  manifolds  also  work  only  for  a  cross¬ 
type  gas-flow  configuration.  Internal  manifolds  are  often 
advantageous  not  only  because  of  better  sealing,  but  also 
because  the  stack  is  not  sensitive  to  a  change  in  height  due  to 
shrinking.  In  a  molten  carbonate  fuel-cell  (MCFC)  stack,  for 
example  stack  height  changes  significantly  as  the  temperature 
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rises  from  room  temperature  to  its  operating  condition 
(650  °C)  because  of  shrinking  of  matrix  layers.  In  this  respect, 
internal  manifolds  are  more  practical  than  external  designs. 
Shrinking  may  not  be  significant  for  solid  oxide  or  polymer 
electrolyte  membrane  fuel-cell  stacks,  but  the  internal  mani¬ 
fold  design  is  usually  favored  for  these  fuel-cells  because  it 
allows  more  versatility  in  gas-flow  configuration. 

The  geometrical  design  factors  of  internal  manifold  stacks 
include  manifold  structure,  size,  number  of  manifolds,  over¬ 
all  gas-flow  pattern,  gas  channel  depth,  active  area  for 
electrode  reactions,  and  so  on.  The  shapes  of  manifold  holes 
can  be  circular,  rectangular,  or  oblong.  The  cross-sectional 
area  of  the  holes  is  important  because  it  determines  the  linear 
velocity  of  gas-flow  through  the  manifolds  for  a  given  flow 
rate  of  the  inlet  gas.  The  overall  gas-flow  pattern  consists  of 
several  different  flow  components,  namely  inlet  manifold 
flow,  flow  through  channels  on  the  electrodes,  and  outlet 
manifold  flow.  Usually,  the  flow  pattern  can  be  either  a 
U-shape  (reverse  flow)  where  the  outlet  gas  flows  in  opposite 
direction  to  the  inlet  gas,  or  a  Z-shape  (parallel  flow)  where 
the  directions  of  the  inlet  and  outlet  gas-flows  are  the  same 
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Nomenclature 

Ax  cross-section  area  of  gas-flow  channel  on 

electrodes  (m2) 

D  diameter  of  pipe  (m) 

Dh  hydraulic  diameter  of  manifold  hole  (m) 

/  fanning  friction  factor 

Ft  force  exerted  by  fluid  during  motion  (N  m  2) 

g  gravitational  acceleration  constant 

(=9.8  m  s“2) 
h  height  (m) 

k  permeability  (m2) 

Kf  geometrical  pressure-loss  coefficient 

/wf  friction  loss  in  mechanical  energy  balance 

(m2  s'2) 

L  length  of  pipe  (m) 

Le  equivalent  cell-to-cell  length  in  manifold  flow 

region  (m) 

Lx  length  of  electrodes  in  fuel-cell  stack  (m) 

N  number  of  cells  in  a  stack 

Pexit  exit  pressure  at  outlet  manifold  (N  m-2) 

AP(i )  pressure-change  across  the  ith  cell  (N  m-2) 
R(i)  local  gas  dynamic  pressure  at  the  ith  cell 
(N  m“2) 

qx{i )  gas-flow  rate  in  the  ith  cell  channel  (m3  s_1) 

QF  total  feed  gas-flow  rate  (m3  s_1) 

Q(i)  cumulative  gas-flow  rate  in  the  ith  manifolds 

(m3  s-1) 

Red)  Reynolds  number  at  the  ith  cell 

w(i)  gas  linear  velocity  at  the  ith  cell  (m  s-1) 

Greek  letters 

ad)  pressure  distribution  factor  for  ith  cell,  defined 
in  Eq.  (21) 

e  tolerance  for  calculation  error 

i]  viscosity  (kg  m~'  s-1) 

p  density  (kg  m-3) 

ts  shear  stress  (N  m~2) 

Superscripts 

adj  adjusted  value  during  inner  loop  of  calculation 
new  newly  calculated  value  for  next  iteration  step 

Subscripts 

in  inlet  manifolds 

out  outlet  manifolds 

ref  reference  cell 

x  coordinate  of  flow  through  cell  channel 

y  coordinate  of  flow  through  manifolds 


(Fig.  1).  In  both  cases,  the  inlet  and  outlet  manifold  flows  are 
perpendicular  to  the  gas-flow  through  the  channels  on  the 
electrodes.  Other  types  of  flow  pattern  may  also  be  possible. 

ft  is  often  emphasized  in  fuel-cell  stack  design  that 
pressure -change  in  manifolds  should  be  much  lower  than 


Gas  out 


Parallel  flow 
(Z-shape) 


Gas  in 


Fig.  1.  Schematic  of  two  different  types  of  gas-fiow  pattern  in  fuel-cell 
stacks. 


that  in  the  gas  channels  on  the  electrodes  in  order  to  ensure  a 
uniform  flow  distribution  among  cells  piled  in  a  stack.  This 
rule-of-thumb  is  often  satisfied  when  the  number  of  cells  in  a 
single  stack  is  not  too  large.  There  is  no  general  design  rule, 
however,  as  to  how  many  cells  are  allowed  to  make  a  single 
stack  with  a  reasonably  good  flow  distribution.  There  have 
not  been  many  large-scale  stacks  tested  with  near  100-cells 
in  a  single  stack.  One  of  the  reasons  for  this  paucity  of  test 
results  is  that  fuel-cell  stack  construction  is  too  costly  to 
undertake  mainly  for  the  purpose  of  a  pilot-scale  engineer¬ 
ing  test.  Theoretical  approaches  are  therefore  favored  in  the 
engineering  design  of  fuel-cell  stacks.  It  is  the  purpose  of 
this  paper  to  suggest  a  method  to  calculate  the  pressure 
variation  for  each  directional  component  of  gas-flow  in 
planar  type  fuel-cell  stacks  in  relation  to  flow  distribution 
among  the  individual  cells.  There  have  been  a  few  related 
studies  in  recent  years  [2-9],  but  we  propose  a  more 
systematic  algorithm  for  calculation.  A  specific  case  is 
illustrated  by  applying  the  proposed  method  for  a  pilot-scale 
fuel-cell  stack. 


2.  Mathematical  model 

2.1.  Number  of  variables  and  degree  of 
freedom  for  analysis 

Two  overall  flow  patterns  in  fuel-cell  stacks  with  planar 
cells  are  depicted  in  Fig.  1.  The  U-shape  gas-flow  may  be 
called  a  ‘reverse  flow’  as  the  outlet  gas-flows  in  opposite 
direction  to  the  inlet  flow.  The  Z-shape  gas-flow  is  a  ‘parallel 
flow’  since  both  inlet  and  outlet  gases  flow  in  the  same 
direction.  Let  us  consider  a  stack  where  cells  are  numbered 
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where  QP  =  total  inlet  feed  flow  rate. 

A/>in(t)  =  -Pin(t')  -  P-m{i  -  1),  (*  =  2,3,..., TV)  (4) 

In  outlet  manifolds  of  a  U-shape  gas-flow  stack,  Qout(i) 
and  Al\mt(i)  are  defined  as  follows: 

i 

Gout(0  =  ^  'fls(j)  —  Cout(i  1)  fr(J)  (5) 

j=  1 

APout(i)  =  PoutO'  -  1)  -  ^out(i);  (*  =  2, 3, . . .  ,N) 

(6) 

In  outlet  manifolds  of  a  Z-shape  gas-flow  stack,  the  same 
variables  are  given  by: 

N 

Gout(i)  =  =  fiout(*  +  !)  +  #*(0  (7) 

7=1 
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Fig.  2.  Designation  of  pressure  and  flow  variables  for  mathematical 
model:  (a)  reverse  flow,  (b)  parallel  flow. 


from  the  highest  (1st)  to  the  lowest  (/Vth)  position  and  the 
feed  gas  enters  the  Mh  cell  as  shown  in  Fig.  2a  for  a  reverse 
flow  and  in  Fig.  2b  for  a  parallel  flow.  The  gas-flow  in  each 
cell  through  the  channel  area  is  denoted  by  qx(i).  The 
dynamic  pressure  drop  from  the  inlet  to  the  outlet  of  a 
gas  channel  is  denoted  by  A Px(i)  which  must  be  a  function  of 
cell  flow  rate,  i.e. 

A Px{i)  =  Pm(i)  -  POut(0,  («'  =  1, 2, 3, . . . , N)  (1) 

In  inlet  manifolds,  the  cumulative  flow  rate,  and  the 

local  pressure  drop,  APin(i),  are  defined  as  follows: 

i 

Gin(i)  =  L ^(7 )  =  Gin(*  -  !)  +  £c(0  (2) 

j=  1 

Gin  (AO  =  Qf  (3) 


AFout(0  =  P<mt(i)  —  P<mt(i  —  1),  (i  =  2,3, . . .  ,N)  (8) 

From  Fig.  2,  note  that  the  total  inlet  flow  rate  is  Qm(N )  in 
both  flow  types  while  the  total  outlet  flow  rate  is  Qout(N )  in  a 
U-shape  flow  and  Gout(l)  in  a  Z-shape  flow. 

For  an  A'-cell  stack,  there  are  5 N  +  1  variables  involved  to 
represent  local  pressure  and  flow  rate,  viz.  Pm(i),  P(mt(i), 
Gin(0,  Gout(f),  q*(i),  and  QF.  Among  them,  Q  Ji)  and  Gout(0 
are  related  directly  to  qx(i)  from  the  2N  mass  balance 
equations  (Eqs.  (2),  (5)  and  (7)).  A  feed  flow  rate  ( QP ) 
should  be  given  as  one  of  the  operating  parameters.  To  solve 
the  given  set  of  equations,  therefore,  requires  3 N  indepen¬ 
dent  equations  that  relate  local  pressure  and  cell  flow  rate, 
qx( i).  Two  relations  are  pre-determined  because  Qm (N)  is 
equal  to  the  feed  flow  rate  (QP)  and  the  exit  pressure  must  be 
known.  The  exit  pressure  is  Pout(N )  for  a  U-shape  flow  and 
Poutd)  for  a  Z-shape  flow.  The  exit  pressure  is  usually  set  to 
zero  for  dynamic  pressure  calculations.  There  should  be 
(N  —  1 )  equations  for  APin(i)  and  another  (N  —  1)  equations 
for  APout(i).  In  addition,  A  equations  can  be  obtained  from  a 
relationship  between  cell  flow  rate,  qji),  and  pressure  drop, 
A Px(i).  Once  all  these  relationships  are  defined,  the  degree  of 
freedom  becomes  zero  and  thus  the  system  can  be  solved  to 
obtain  flow  and  pressure  distributions  inside  the  manifolds 
and  the  cell  channels. 

2.2.  Relation  for  cell  flow  rate  and  pressure  drop 

The  gas-flow  in  cell  channels  can  be  approximated  to  a 
flow  through  a  conduit  for  which  many  fluid  dynamics 
relations  are  available  as  a  function  of  the  Reynolds  number. 
Therefore,  a  flow  model  can  be  formulated  to  describe  the 
gas-flow  in  cell  channels.  During  operation  under  a  current 
load,  the  mass  flow  rate  varies  from  channel  inlet  to  channel 
outlet  as  a  result  of  consumption  or  generation  by  cell 
reactions.  Such  changes  in  mass  flow  should  be  accommo¬ 
dated  in  the  flow  model  when  gas  utilization  is  high.  When 
gas  utilization  is  low,  however,  such  a  change  of  mass  is  not 
significant.  The  change  of  mass  is  not  included  in  the  model 
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explained  below,  but  it  can  be  easily  added  to  the  model 
equations  with  some  modification.  The  treatment  of  mass 
change  due  to  cell  reactions  has  been  reported  in  a  study  of 
fuel-cell  performance  models  [10]. 

The  gas-flow  area  on  electrodes  is  often  not  straight 
channels  but  includes  some  regularly  shaped  obstacles  for 
effective  flow  dispersion.  These  obstacles  cause  more  flow 
resistance  and,  consequently,  a  higher  pressure  drop  than 
straight  open  channels.  The  pressure  drop  in  the  gas  channels 
is  therefore  strongly  dependent  on  the  shape  and  the  config¬ 
uration  of  the  flow  obstacles.  A  system-dependent  parameter 
is  necessary  for  the  flow  model  in  cell  channels  to  absorb 
the  effect  of  flow  obstacles  in  gas  channel  area.  Since  the  flow 
in  a  cell  channel  is  usually  in  a  laminar  region,  a  linear 
relation  of  A P  with  flow  rate  should  be  reasonably  satis¬ 
factory.  A  model  of  the  cell  channel  flow  can  be  derived 
from  the  Hagen-Poiseullie  equation,  Darcy’s  law,  or  others. 
Here,  Darcy’s  law  is  used  in  which  permeability  ( k )  is  the 
parameter  that  accounts  for  the  flow  resistance  of  the  channel 
area. 


d  P 
dx 


(9) 


Assuming  gas  viscosity  and  velocity  are  constants,  the 
pressure  drop  from  the  inlet  to  the  outlet  is  represented  as 
follows: 


A  Px(i) 


P  in  ( i)  Hout(i) 


(10) 

(ID 


It  is  necessary  to  determine  the  permeability  from  mea¬ 
surement,  and  then  the  pressure  drop  from  each  cell  channel 
is  obtained  directly  from  the  linear  gas  velocity  at  the  ;'th 
cell,  ux{i). 


2.3.  Relation  for  manifold  pressure  and  flow  rate 


The  flow  in  manifolds  is  described  quite  differently  from 
the  flow  in  cell  channels.  In  inlet  manifolds,  the  gas-flow  rate 
decreases  from  the  feeding  position  to  the  other  end.  In  outlet 
manifolds,  the  gas  flow  rate  increases  as  the  flow  from  each 
cell  channel  is  sequentially  added  to  the  main  flow.  The 
overall  gas-flow  in  manifolds  is  characterized  by  manifold 
axial  flow  (uy)  in  combination  with  local  dynamic  pressure 
(P).  A  mathematical  model  for  the  gas  pressure  in  manifolds  is 
developed  from  an  overall  mechanical  energy  balance,  or  the 
engineering  Bernoullie  equation.  Local  velocities  in  the 
mechanical  balance  equation  are  related  to  mass  or  volumetric 
flow  rates,  which  must  be  conserved  from  mass  balances,  i.e. 


A  Uy2 

~T~ 


+  gAh  + 


A  P 
P 


(\V'I 


(12) 


The  friction  loss  term  (Zwf)  includes  wall  friction  loss  with  a 
wall  friction  coefficient  (f)  which  is  defined  as  a  ratio  of 


shear  stress  to  inertial  stress,  i.e. 


shear  stress 
inertia  stress 


(1/2  )puy2 


(13) 


The  force  exerted  by  fluid  during  a  motion  in  a  pipe  is 
expressed  by  shear  stress  and  dynamic  pressure-loss  by 
definition,  i.e. 


Fr 


=  xsnDL  =  A  P 


(14) 


By  combining  Eqs.  (13)  and  (14),  the  wall  friction  loss  is 
expressed  in  the  most  frequently  available  form  as  follows: 


A-Ploss  -  2 pUy  f  I 


D 


(15) 


Additionally,  there  is  a  friction  loss  due  to  geometrical 
effects  such  as  fittings  or  valves.  In  the  gas  manifolds  of  fuel¬ 
cell  stacks,  there  are  some  geometrical  effects  not  from 
fittings  or  valves  but  from  the  flange-shape  structure  at  the 
junction  of  the  manifold  axial  flow  and  each  cell  channel 
flow.  A  schematic  illustration  of  the  structure  used  in  the 
stack  considered  in  this  study  is  given  in  Fig.  3. 

Such  a  geometrical  effect  is  generally  estimated  from  the 
inertia  force  multiplied  by  an  empirical  constant  (Kf)  which 


Fig.  3.  Shape  of  manifold  holes  in  bipolar  plates  of  a  molten  carbonate 
fuel-cell  stack. 
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is  known  for  various  types  of  fittings  or  valves.  The  com¬ 
bined  friction  loss  (wall  +  geometry)  is  then  obtained  as  a 
function  of  linear  velocity,  namely: 


/wf  =  wall  friction  loss  +  fitting  loss 


uv 

(16) 


The  overall  mechanical  balance  is  now  expressed  with  an 
equivalent  length  of  manifold  holes  from  cell  to  cell  ( Le )  and 
a  hydraulic  diameter  (Z)H),  i.e. 


A  Uy2 


+  gAh  + 


A  P 
P 


(17) 


The  potential  energy  change  term  (gAh)  can  be  ignored  as 
it  is  usually  smaller  than  the  others  by  an  order  of  magnitude 
at  least.  In  inlet  gas  manifolds,  the  kinetic  energy  term 
(Ally2)  contributes  a  negative  value  of  pressure  drop  which 
indicates  pressure  recovery.  This  occurs  when  the  mass  flow 
in  inlet  manifolds  is  reduced  from  the  feeding  plate  to  the 
end  plate,  which  result  in  an  increase  in  linear  velocity. 
When  the  friction  loss  is  relatively  small  compared  with  the 
pressure  recovery,  the  overall  A P  is  negative  and  pressure 
increases  along  the  main  flow  direction  of  the  inlet  mani¬ 
folds.  The  reverse  is  true  when  there  is  a  relatively  large 
friction  loss. 

The  equation  at  each  junction  of  the  manifolds  is  shown 
below.  This  equation  can  apply  for  both  inlet  and  outlet 
manifolds  by  replacing  the  pressure  term  (P)  by  Pm  or  Pout, 
i.e. 


A P® 
P 


P(i)-P(i-1) 

P 

Uy(i  ~  l)2  —  Uy(l)~ 

2 


+  2 fuy(i)1 


+  Kf 


Uy(i)2 

2 


(18) 


While  the  gas-flow  in  the  cell  channel  is  mostly  in  the 
laminar  region,  gas-flow  in  manifolds  can  be  either  laminar 
or  turbulent.  The  coefficient  of  wall  friction  loss  can  be 
obtained  from  a  known  correlation  [11]: 

/  =  ,  Re  <  2100  (Laminar  flow)  ( 1 9) 

/  =  0.079  Re~0-25, 

Re  >  2100  (Blasius  equation  for  turbulent  flow)  (20) 


2.4.  Solution  algorithm 

The  calculation  starts  with  an  initially  guessed  set  of  cell 
flow  rates,  qx(i).  The  easiest  way  to  start  is  with  an  equally- 
divided  cell  flow  rate,  i.e.  the  feed  flow  rate  divided  by  the 
number  of  cells.  With  a  given  set  of  cell  flow  rates,  the 
cumulative  local  flow  rates  in  manifolds,  Qin(i)  and  Qoul(i), 
are  obtained  from  Eqs.  (2),  (5)  and  (7).  Local  velocities  in 


manifolds  are  obtained  from  the  flow  rates  and  these  are  then 
used  to  calculate  the  pressure-change  in  the  outlet  manifolds 
from  Eq.  (18)  and  the  local  dynamic  pressure  in  the  outlet 
manifolds  (Pout(0)  from  Eqs.  (6)  and  (8).  The  next  step  is 
to  calculate  Pjn(i)  from  the  (N  —  1)  relations  for  Pin(i)  and 
uy(i)  (Eq.  (18)).  Since  there  are  only  (N  —  1)  equations  for 
N  pressure  variables,  it  is  necessary  to  specify  one  local 
pressure  value  at  the  inlet  manifolds.  This  is  called  the 
reference  inlet  pressure,  Pin(ref),  and  is  easily  obtained  from 
the  pressure  drop  in  a  reference  cell,  APv(ref).  The  reference 
cell  can  be  chosen  arbitrarily.  Here,  No.  1  cell  (top  of  the 
stack)  is  taken  as  the  reference. 

Examination  of  the  variables  and  equations  reveals  that 
there  are  more  equations  for  pressure  than  for  flow  rate.  It  is 
therefore  useful  to  define  a  parameter  which  utilizes  a 
redundant  set  of  pressure  equations.  A  basic  solution  strategy 
is  then  to  use  this  parameter  for  adjustment  of  cell  flow  rates 
and  continue  until  convergence  reaches.  Here,  the  redundant 
parameter  is  defined  as: 


_  Pin(i)  ~  Pout(0 
A/5*  (ref) 


(21) 


where  APx(ref)  is  the  pressure  drop  in  the  reference  cell. 
When  all  the  flows  in  cells  are  equal,  the  ct(i)  values  are  all 
unity.  A  large  deviation  of  a(i)  values  from  cell-to-cell 
means  the  total  gas  flow  is  not  uniformly  distributed  to 
cells.  Therefore  a.(i)  may  be  called  a  ‘flow  distribution 
factor’ . 

Some  adjustment  correlation  is  needed  to  make  a  better 
guess  for  the  next  calculation  in  each  iteration  step.  The 
adjustment  correlation  is  developed  from  the  definition  of 
the  flow  distribution  factor,  namely: 


—  AP  (ref)  'WOl  @2) 

1—1  1  1=1 

By  substituting  the  pressure  drop  terms,  Pm(i)  -  Pout®  in 
Eq.  (22)  by  the  equation  for  cell  flow  rate  (Eq.  (11))  the 
following  formula  is  obtained  for  a  new  value  of  the  pressure 
drop  in  the  reference  cell: 


AP"ew(ref) 


Qf  1  i]Lx 

EL«(0^ 


(23) 


With  this  new  value  of  pressure  drop  in  the  reference  cell, 
the  inlet  manifold  pressure  values  are  all  adjusted.  The  flow 
distribution  factor  should  be  adjusted  as  well.  New  cell  flow 
rates  are  then  calculated  from  the  adjusted  inlet  manifold 
pressure  values.  Note  that,  the  outlet  manifold  pressure  is 
not  adjusted  because  the  exit  pressure  is  known  and  is 
equivalent  to  the  outlet  reference  pressure.  Once  new  cell 
flow  rates  are  determined,  they  are  compared  with  the 
initial  guess  or  the  previously  calculated  cell  flow  rates. 
The  calculation  continues  with  a  new  set  of  cell  flow  rates 
until  the  qx(i)  values  are  close  to  the  old  values  within  an 
error  limit.  The  entire  solution  algorithm  is  summarized  in 
Appendix  A. 
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3.  Experimental 

The  gas  pressure  from  a  10-cell  MCFC  stack  was  mea¬ 
sured.  The  cell  size  was  55  cm  length  and  110  cm  width 
which  resulted  in  an  active  cell  area  of  6050  cm2.  The  cell 
component  materials  and  other  specifications  have  been 
discussed  elsewhere  [12].  The  fuel  and  oxidant  gases  were 
fed  through  internal  gas  manifolds.  Both  the  cathode  and 
anode  gas  streams  flow  in  the  same  direction  from  the  inlet  to 
the  outlet  (co-flow).  The  overall  gas-flow  pattern  is  a  reverse 
type  (U-shape).  The  gas  pressure  was  measured  from  inlet 
and  outlet  headers  of  the  stack  at  different  gas-flow  rates. 
The  overall  plane  view  of  separator  plates  and  a  schematic 
view  of  the  manifold  cross-section  are  shown  in  Fig.  3. 

4.  Results  and  discussion 

4.1.  Pressure  drop  across  gas  channels  in  cell  area 

The  measured  gas  pressure  drop  is  plotted  as  a  function  of 
flow  rate  in  Fig.  4.  The  cathode  and  anode  plates  have  the 
same  geometrical  structure  with  the  same  dimensions.  It  is 
therefore  reasonable  that  the  pressure-drop  data  for  the 
anode  and  the  cathode  in  Fig.  4  are  almost  identical.  Linear 
fitting  of  data  results  in  two  close  straight  lines  for  the  anode 
and  the  cathode,  respectively.  The  small  difference  between 
the  anode  and  cathode  lines  is  probably  a  result  of  other 
effects  such  as  different  gas  viscosity  and  different  electrode 
materials.  The  linear  behavior  of  pressure  drop  indicates  that 
the  gas-flow  is  in  the  laminar  region.  It  should  be  noted  that 
the  data  was  obtained  from  a  10-cell  stack.  Though  the  data 
were  measured  between  inlet  and  outlet  headers  that  pre¬ 
serve  the  effect  of  manifold  pressure,  the  effect  of  gas 
pressure-change  in  manifolds  should  not  be  significant  in 


Fig.  4.  Pressure  drop  across  cell  channels  measured  from  55  cm  length 
MCFC  cells. 


this  short  stack  of  only  10-cells.  The  pressure  drop  in  Fig.  4 
is  caused  mostly  by  the  flow  resistance  in  cell  channels,  A Px. 
From  the  relation  based  on  Darcy’s  law,  the  permeability  can 
be  used  to  characterize  the  flow  resistance  in  cell  channels. 
The  permeability  for  anode  and  cathode  cell  channels  was 
calculated  from  the  data,  and  the  results  are  listed  in  Table  1 . 
The  values  for  the  anode  and  cathode  are  very  close,  as 
expected. 

When  a  stack  is  composed  of  much  more  than  10-cells, 
the  pressure-change  in  the  manifolds  may  become  signifi¬ 
cant  given  the  size  and  shape  of  the  manifold.  While  a  short 
stack  is  easy  to  set  up  and  test,  a  stack  with  a  large  number  of 
cells  is  always  difficult  to  test.  The  pressure  drop  from  the 
cell  channels  can  be  measured  even  with  a  single  cell,  and 
therefore  it  is  relatively  easy  to  determine  the  flow  resistance 
of  the  cell  channels,  as  undertaken  in  this  study.  The 
pressure-change  in  the  manifolds  of  a  large  stack  is  possible, 
but  costly.  In  the  next  section,  the  manifold  flow  is  described 
with  the  flow  model  (Eq.  (18))  and  local  pressure  values  are 
predicted.  A  stack  model  with  100-cells,  each  with  an 
identical  size  of  active  area  to  the  10-cell  stack  described 
above,  is  formulated.  It  should  be  considered  a  scale-up  from 
10  to  100-cells.  Since  the  active  cell  area  is  the  same,  the 
flow  resistance  from  the  10-cell  stack  data  applies  for  cell 
channels  of  a  100-cell  stack.  Thus,  it  is  possible  to  inves¬ 
tigate  the  change  of  pressure  in  manifolds  and  its  effect  on 
flow  distribution.  Both  reverse  and  parallel  flow  types  are 
considered.  Cells  are  numbered  from  the  highest  to  the 
lowest  (see  Fig.  2).  The  feed  gas  enters  the  lowest  plate 
(Ath),  and  the  exhaust  gas  exits  from  the  /Vth  plate  in  a 
reverse  flow  and  from  the  1st  plate  in  a  parallel  flow,  as 
shown  in  Fig.  2.  The  specification  of  this  stack  and  the 
simulation  conditions  are  summarized  in  Table  1. 

4.2.  Analysis  of  flow  in  manifolds 

The  manifold  flow  model  developed  in  this  study 
(Eq.  (18))  is  a  relation  between  manifold  gas  pressure  iP) 
and  velocity  ( uy )  with  irreversible  pressure-loss  character¬ 
ized  by  a  wall  friction  loss  coefficient  (f)  and  a  geometrical 
loss  coefficient  (Kf).  Besides  the  two  pressure-loss  terms,  the 
kinetic  energy  change  term  also  influences  pressure-change 
in  manifolds  because  mass  change  occurs  at  the  junction  of 
the  manifold  and  the  cell  channels.  The  effect  of  such  a 
kinetic  energy  change  is  pressure  recovery  in  inlet  manifolds 
and  pressure -loss  in  outlet  manifolds.  Therefore,  there  are 
three  different  contributions  to  pressure-change  in  mani¬ 
folds,  as  noted  in  Eq.  (24).  The  wall  friction  loss  is  a  pressure 
drop  caused  by  viscous  stress  to  the  smooth  walls  of  the  flow 
length.  The  pressure  recovery  is  related  with  mass  flow 
change  and,  consequently,  velocity  change.  The  geometrical 
loss  exists  because  manifold  holes  are  not  exactly  straight 
pipes  with  smooth  walls.  Thus, 

Apnanifolds  —  APmass  change  T  A/J|'|-jctjon  ]oss 

=  AP  mass  +  A  P  wan  +  A  P  ge0  (24) 
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Table  1 


Specification  of  stack  and  its  standard  operating  conditions 


Variable 

Value 

Number  of  cells 

100 

Gas-flow  type 

Reverse  (U)  or  parallel  (Z) 

System  total  pressure  (atm) 

3 

Stack  heating  temperature  (°C) 

650 

Electrode  size 

0.55  m  X  1.10  m 

Effective  gas  channel  depth  (m) 

0.0012 

Cell-to-cell  distance  in  manifold  axis  (m) 

0.005 

Base  load  current  density  (mA  cm-2) 

125 

Cathode  gas  O2  utilization 

0.3 

Anode  gas  H2  utilization 

0.8 

Cathode  inlet  gas  composition 

02/N2/C02  =  15/55/30 

Anode  inlet  gas  composition 

H2/C02/H20  =  72/18/10 

Cathode  molar  feed  flow  rate  per  cell  (mole  s-1) 

0.0432 

Anode  molar  feed  flow  rate  per  cell  (mole  s_  ) 

0.00675 

Cathode  gas  viscosity  (kg  m_I  s-1) 

3.6  x  10~5 

Anode  gas  viscosity  (kg  m-1  s_  ) 

2.9  x  l(r5 

Cathode  gas  density  (kg  m-3) 

1.41 

Anode  gas  density  (kg  m_  ) 

0.54 

Cathode  gas  permeability  in  cell  channel  (m2) 

1.863  x  10~8a 

Anode  gas  permeability  in  cell  channel  (m2) 

1.616  x  10~8a 

Number  of  cathode  gas  manifold  holes 

5  holes  (inlet)  — > 

4  holes  (outlet) 

Number  of  anode  gas  manifold  holes 

4  holes  (inlet)  — ► 

5  holes  (outlet) 

Hydraulic  diameter  of  cathode  manifold  holes 

0.0786  m  (inlet)  - 

->  0.0691  m  (outlet) 

Hydraulic  diameter  of  anode  manifold  holes 

0.0691  m  (inlet)  - 

->  0.0786  m  (outlet) 

a  Obtained  based  on  data  of  Fig.  4. 


The  pressure-change  due  to  the  mass  change  ( A/Jnlass)  is  a 
reversible  energy  change  caused  by  different  gas  velocities 
before  and  after  the  junction  of  manifolds.  The  pressure-loss 
due  to  wall  friction  (A/\van)  is  estimated  with  a  wall  friction 
loss  coefficient  from  correlation  (Eqs.  (19)  and  (20)).  The 
pressure-loss  due  to  geometrical  effect  (A Pgeo)  is  estimated 
with  a  geometrical  coefficient.  Though  some  geometrical 
pressure-loss  coefficients  are  available  for  various  fittings  or 
valves,  none  fit  a  specific  shape  of  gas-flow  manifolds  such 
as  that  shown  in  Fig.  3.  The  shape  of  the  manifold  junction 
can  be  approximated  to  some  known  types  of  fitting  or  valve. 
In  Table  2,  several  different  cases  of  approximation  are 
suggested  and  Kf  values  are  compared  [13].  There  are  also 
different  approaches  to  calculate  pressure  and  flow  distribu¬ 
tion  for  other  types  of  application  of  manifolds  [14]. 

A  typical  plot  of  Reynolds  number  versus  cell  number  is 
shown  in  Fig.  5.  The  y-coordinate  is  plotted  in  a  logarithmic 

Table  2 


Geometrical  pressure-loss  coefficients  from  different  approximations  [13] 


Geometry 

Kf  formula 

Kf  value3 

Straight  pipe 

0.0 

0.0 

Rounded  entrance 

0.05 

0.05 

Sudden  contraction 

0.45  (1  -  f) 

0.2 

Fully-open  gate-valve 

0.2 

0.2 

Orifice  (sharp-edge) 

2.7  (1  -/?)(!-  p2)  a/p2) 

2.679 

Here  P  =  (smaller  cross-sectional  area)/(larger  cross-sectional  area). 
3  For  fuel-cell  stack  in  this  study,  p  =  0.556. 


scale.  When  the  plot  is  made  in  a  linear  scale,  it  is  almost 
linear  with  the  variation  of  Re  in  the  cells.  Though  Fig.  5 
shows  the  case  of  a  reverse  flow  only,  the  Reynolds  number 
in  a  parallel  flow  is  qualitatively  the  same  for  Rem  and  Rex, 
while  R<?0ut  should  be  symmetrical  with  respect  to  the  center- 
height  cell  plane  (50th).  The  range  of  Reynolds  numbers 
for  cell  channel  flow  ( Rex )  is  below  the  order  of  102  in  Fig.  5 
and  is  certainly  in  the  laminar  flow  region.  The  range  of 
Reynolds  numbers  for  manifold  flow  is  placed  in  both  the 
laminar  and  turbulent  regions.  The  manifold  flow  curves 
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cross  the  line  of  criterion  for  flow  type  (Re  =  2100)  at  a  cell 
number  of  about  10-20.  This  means  that  turbulent  flow  is 
expected  if  the  number  of  cells  exceed  20  in  this  stack 
configuration  and  manifold  design.  For  laminar  flow,  the 
pressure  drop  is  linearly  proportional  to  the  gas  velocity 
and  this  usually  creates  a  less  significant  pressure-change 
compared  with  turbulent  flow  where  the  pressure  drop  is 
amplified  with  increasing  velocity,  for  example  (velocity)1'75 
from  the  Blasius  equation.  For  a  100-cell  stack,  much  of  the 
manifold  flow  is  turbulent  and  this  will  create  a  larger 
pressure  drop  than  for  laminar  flow. 

4.3.  Pressure-loss  and  recovery  in  manifolds 

With  the  Kf  values  in  Table  2,  the  flow  model  is  solved 
according  to  the  algorithm  shown  in  Appendix  A.  The  three 
pressure-change  terms  in  Eq.  (24)  are  plotted  for  inlet 
manifold  as  a  function  of  cell  number  in  Fig.  6.  A  positive 


A P  indicates  pressure-loss,  whereas  a  negative  A P  indicates 
pressure  recovery.  When  there  is  no  geometrical  effect 
(Kf  =  0),  the  pressure-loss  is  negligible  compared  with 
the  pressure  recovery  caused  by  mass  decrease  in  inlet 
manifolds  (Fig.  6a).  As  Kf  increases  (Fig.  6a-d),  the  magni¬ 
tude  of  the  geometrical  pressure-loss  increases  dramatically 
while  the  magnitude  of  the  pressure  recovery  shows  only  a 
small  increase.  Pressure-loss  by  wall  friction  is  always  much 
smaller  than  both  the  other  effects  regardless  of  the  value  of 
Kf.  The  two  dominant  contributions  of  pressure-change  in 
manifolds  are  obviously  from  mass  change  and  geometry, 
and  their  relative  influence  depends  greatly  on  geometrical 
shape  that  determines  a  Kj  value. 

The  data  in  Fig.  6  demonstrate  how  manifold  geometry 
affects  the  flow  in  manifolds.  The  loss  coefficient  of  0.2 
is  for  a  fully-open  gate-valve  that  has  a  relatively  smaller 
flow  resistance  compared  with  most  other  fittings  or  valves. 
Despite  such  a  relatively  small  geometrical  flow  resistance, 


Fig.  6.  Analysis  of  pressure-loss  and  recovery  in  inlet  manifolds  at  different  values  of  geometrical  pressure-loss  coefficient:  (a)  Kf  —  0.0;  (b)  Kf  —  0.05; 
(c)  Kf  =  0.2;  (d)  Kf  =  2.679. 
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its  effect  on  the  distribution  of  manifold  pressure  is  sig¬ 
nificant,  as  shown  in  Fig.  6.  This  explains  why  the  manifold 
is  a  crucial  part  in  the  design  of  large-scale  fuel-cell  stacks. 

It  is  most  difficult  to  determine  the  Kj  value,  since  fuel¬ 
cell  stacks  can  be  designed  in  many  different  ways.  Because 
flow  resistance  must  be  affected  by  both  the  structure  and  the 
size  of  the  stack,  it  is  safe  to  assume  that  every  stack  has  a 
different  flow  resistance  in  the  manifolds  and,  therefore,  a 
different  Kf  value.  If  we  take  an  example  from  Fig.  3,  there 
are  no  smooth  walls  in  the  gas-flow  through  manifolds. 
These  manifold  holes  are  more  like  pipes  with  rough 
surfaces  on  the  walls  and  with  a  series  of  many  orifices 
in  the  flow  direction.  Though  some  pressure-loss  coefficients 
are  available  in  literature,  most  of  them  are  for  flow  in  pipes 
with  smooth  walls  and  for  very  well-described  geometry. 
The  manifold  shape  can  still  be  assumed  to  be  one  of  the 
known  geometry  types  of  fittings  or  valves,  but  some  pre¬ 
cision  of  the  model  may  have  to  be  sacrificed  when  such 
approximations  are  made. 


4.4.  Overall  pressure  variation  in  manifolds 

Combining  the  pressure  drop  in  cell  channels  and  each 
contribution  of  pressure-change  in  inlet  and  outlet  mani¬ 
folds,  the  overall  pressure  distribution  profiles  along  the 
gas-flow  path  in  the  stack  is  obtained,  shown  in  Fig.  7.  For 
comparison,  the  plots  are  made  for  the  cases  of  Kf  =  0.0  and 
0.2.  Both  reverse  and  parallel  flows  are  considered  for  each 
Kf  value.  When  there  is  no  geometrical  pressure-loss 
(Kf  =  0.0),  there  is  virtually  no  noticeable  change  in  pres¬ 
sure  from  cell-to-cell,  i.e.  the  pressure  curves  are  almost  flat 
(Fig.  7a  and  b).  More  careful  examination,  however,  reveals 
that  the  inlet  and  outlet  pressure  curves  are  both  increasing  in 
parallel  in  a  reverse  flow  (Fig.  7a)  and  diverge  from  the 
lowest  cell  No.  (100)  to  the  highest  cell  No.  (1)  in  a  parallel 
flow  (Fig.  7b).  When  there  is  some  contribution  from  the 
geometrical  effect  (Kf  =  0.2),  the  inlet  and  outlet  pressure 
curves  more  clearly  converge  from  the  lowest  to  the  highest 
cell  in  a  reverse  flow  (Fig.  7c),  and  decrease  in  parallel  in  a 
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Fig.  7.  Overall  pressure  distribution  in  100-cell  fuel-cell  stack:  (a)  reverse  flow  with  Kf  —  0.0;  (b)  parallel  flow  with  Kf  —  0.0;  (c)  reverse  flow  with 
Kf  —  0.2;  (d)  parallel  flow  with  Kf  —  0.2. 
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parallel  flow  (Fig.  7d).  These  four  different  pressure  curves 
describe  most  general  types  of  pressure  distribution  in  the 
manifolds  of  fuel-cell  stacks.  The  slope  of  each  pressure 
curve  will  change  as  either  the  manifold  size  or  the  manifold 
shape  changes. 

Prediction  of  pressure  curves  is  important,  as  flow  dis¬ 
tribution  into  cells  is  directly  related  to  pressure  variation. 
For  the  100-cell  MCFC  stack  considered  in  this  study, 
lowering  the  pressure -loss  (lower  Kf)  appears  most  bene¬ 
ficial  because  Fig.  7  implies  there  is  almost  even  pressure 
distribution  at  Kf  =  0.0.  This  stack  apparently  has  manifold 
holes  that  are  sufficiently  large  to  prevent  significant  pressure 
recovery  or  drop  from  mass  change.  In  general,  it  is  more 
desirable  to  have  parallel  curves  of  inlet  and  outlet  pressure 
because  the  flow  rate  at  the  /th  cell  ( qx(i ))  is  proportional  to 
the  pressure  difference  between  inlet  and  outlet. 

4.5.  Flow  distribution  in  cells 

From  the  pressure  variation  in  the  manifolds,  the  mass 
flow  to  each  cell  is  calculated.  The  overall  flow  distribution 
with  various  geometrical  loss  coefficients  in  a  reverse  flow 
stack  is  presented  in  Fig.  8.  In  this  plot,  each  cell  flow  rate  is 
divided  by  the  average  flow  rate  over  100-cells.  Even  though 
it  is  not  possible  to  confirm  which  Kf  value  would  best 
approximate  the  manifold  flow  resistance,  at  least  an  indica¬ 
tion  is  given  of  the  flow  distribution  and  the  extent  to  which  it 
can  change  with  increasing  flow  resistance.  The  pressure 
curves  of  Fig.  7  show  that  the  flow  distribution  is  very 
uniform  in  the  complete  absence  of  a  geometrical  effect. 
The  maximum  and  minimum  flow  deviation  is  less  than  4% 
from  the  average,  in  the  case  of  Kf  =  0.  The  flow  distribution 
with  increasing  geometrical  effect  is  still  reasonably  accep¬ 
table  until  Kf  =  0.2.  At  Kf  =  0.2,  the  maximum  flow  devia¬ 
tion  is  about  30%  from  the  average  and  the  minimum  flow 
deviation  is  less  than  10%  from  the  average.  Obviously,  the 


Fig.  8.  Flow  distribution  to  cells  at  different  values  of  geometrical 
pressure-loss  coefficient. 


Fig.  9.  Comparison  of  flow  distribution  to  cells  in  reverse  flow  and 
parallel  flow  for  cathode  and  anode  gases  under  standard  stack  operating 
condition  with  Kf  —  0.2. 

situation  of  Kf  =  2.679  has  to  be  avoided  since  it  results  in  a 
poorly  distributed  gas-flow.  If  some  data  are  available  to 
make  a  close  approximation  of  Kf  then  the  flow  distribution 
can  be  predicted  from  the  data  in  Fig.  8.  If  there  are  no 
experimental  data  for  and  there  is  a  need  to  design  a  large- 
scale  stack,  an  extreme  case  of  non-even  flow  distribution 
may  have  to  be  considered. 

The  overall  flow  distribution  in  a  reverse  flow  and  in  a 
parallel  flow  are  compared  in  Fig.  9  for  an  equal  number 
of  cells  and  sizes  for  Kf  =  0.2.  This  value  of  Kf  was  chosen 
because  it  appears  to  represent  a  fairly  moderate  flow 
resistance  among  the  approximations  listed  in  Table  2. 
Although,  the  two  types  of  flow  display  different  shapes 
of  distribution  curve,  neither  appears  to  be  better  than  the 
other  in  terms  of  flow  uniformity.  In  this  stack,  it  is  clearly 
less  important  whether  the  overall  gas-flow  is  a  U-  or 
Z-shape.  Rather,  a  better  flow  distribution  is  obtained  by 
reducing  the  geometrical  flow  resistance.  This  is  true  only 
for  the  specific  case  selected  in  this  study.  In  some  other 
stack  designs,  it  is  possible  that  either  a  reverse  flow  or  a 
parallel  flow  has  more  benefits  over  other  flow  types  in  terms 
of  better  flow  distribution. 

The  distribution  of  the  cathode  and  anode  gas-flows  in 
the  selected  fuel-cell  stack  are  also  presented  in  Fig.  9.  As 
shown  in  Table  1,  the  anode  gas-flow  rate  is  much  lower  than 
the  cathode  flow  rate  due  to  the  high  utilization  of  fuel. 
Because  of  its  relatively  low  total  flow  rate,  gas  velocities  are 
much  lower  in  anode  gas  manifolds  than  in  cathode  gas 
manifolds.  The  change  of  pressure  in  manifolds  is  therefore 
much  less  significant  at  the  anode  than  at  the  cathode,  and 
this  results  in  an  overall  uniform  flow  distribution  in  the 
anode.  In  general,  the  flow  rate  entering  each  cell  should  be 
larger  than  a  minimally  required  flow  rate.  The  latter  is  the 
rate  of  fuel  supply  to  the  anode  or  oxygen  supply  to  the 
cathode  that  is  completely  utilized  for  a  desired  current  load. 
The  size  of  the  manifolds  should  be  designed  carefully  so 
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that  cell  flow  rates  are  all  above  the  minimum  requirement. 
The  fuel-cell  stack  chosen  for  this  study  was  designed 
satisfactorily  for  its  standard  operation. 


5.  Conclusions 

A  general  model  for  the  calculation  of  pressure  and  flow 
distribution  in  fuel-cell  stacks  has  been  proposed  in  this  study. 
It  is  emphasized  that  pressure  variation  in  fuel-cell  stack 
manifolds  is  influenced  not  only  by  manifold  size  but  also 
by  manifold  geometry.  The  effect  of  manifold  geometry  can 
be  modeled  with  additional  pressure-losses  to  that  experi¬ 
enced  with  a  straight  smooth  wall.  The  pressure-loss  caused 
by  fittings  or  valves  in  the  pipe  flow  can  be  estimated,  but  it  is 
more  desirable  to  determine  a  geometrical  loss  coefficient 
from  experimental  measurements  as  it  must  be  assumed  that 
each  fuel-cell  stack  has  its  unique  design  of  manifold. 

The  manifold  flow  model  combined  with  cell  channel 
flow  has  been  proposed  to  predict  the  pressure  variation  and 
flow  distribution  for  a  100-cell  MCFC  stack.  The  results 
show  that  the  gas-flow  in  cell  channels  is  in  the  laminar 
region  while  the  gas-flow  in  manifolds  is  mostly  turbulent. 
When  the  geometry  of  the  manifold  flow  junction  of  this 
fuel-cell  stack  is  modeled  as  a  fully-open  gate-valve  or 
sudden  contraction  with  a  geometrical  loss  coefficient  of 
0.2,  certain  non-uniformity  of  the  flow  distribution  is 
observed  in  either  a  reverse  flow  (U-shape)  or  a  parallel 
flow  (Z-shape).  The  non-uniform  distribution  is  more  sig¬ 
nificant  at  the  cathode  where  a  large  amount  of  excess  gas 
enters  the  stack,  than  at  the  anode  where  gas  utilization  is 
sufficiently  high  to  allow  a  relatively  small  amount  of  total  gas 
flow.  Overall,  the  model  predicts  that  the  cell  flow  rates  are 
within  a  reasonably  acceptable  range  for  operation  of  the 
chosen  stack  at  its  standard  condition.  It  has  been  demon¬ 
strated  how  the  manifold  design  can  be  evaluated  from 
engineering  models  and  can  be  improved  further  in  this  study. 


Acknowledgements 

This  work  was  supported  by  the  Korea  Electric  Power 
Corporation  (grant  01EC01)  and  the  R&D  Management 
Center  for  Energy  and  Resources  of  the  Korea  Energy 
Management  Corporation. 


Appendix  A.  Algorithm  for  calculation  of  pressure 
and  flow  distributions  in  fuel-cell  stacks 

Step  1 :  Initial  guess  of  cell  flow  rates  and  a  reference  inlet 


pressured  drop 

•J.il)  =  J7:  (i  =  \tON) 

(A.l) 

(ref=1) 

(A. 2) 

Step  2:  Calculate  flow  rates  and  linear  velocity  in  manifolds 
and  Reynolds  numbers 


Gin(0  =  ^ Jlx{j ),  O'  =  1  tO  A) 


j=  1 


GoutO)  =  y ^x(j),  for  reverse  flow 
j=  i 

N 

GoutO)  =  YMj)>  for  parallel  flow 

j=  i 

ux(i),  Uin(i),  uoat(i),  Rex{i),  Rein(i),  Re  out(0 

Step  3:  Calculate  outlet  manifold  pressure 


(A.  3) 


(A. 4) 


^exit  —  0 

(A. 5) 

APout(i)  ^out(^  1)  WOut(0 

P  2 

+  (2/^H+?)"OUt( 

if,  (i  =  2  to  A) 

(A.  6) 

Pout((  —  1)  =  Pout(0  +  APout(O) 

for  reverse  flow 

(A. 7) 

Pout(0  =  Pout{i  ~  1)  +  APout(i), 

for  parallel  flow 

Step  4:  Calculate  inlet  manifold 

pressure 

Pin  (ref)  =  APv(ref)  +  Pout(ref), 

(ref  =  1) 

(A.  8) 

APjnQ)  _  MinQ  ~  l)2  ~  «inQ)2 

P  2 

+  (2/^  +  y)Mm(02,  (i  =  2toA) 

Pin  (0  =  Pin  0  ~  1)  +  A Pin(t) 

Step  5:  Calculate  flow  distribution  factors 

/  .\  ^in(0  PoutW  /•  i  Ar\ 

“(,)=  Af,(,ef)  '  ('=1,0A,) 

Step  6:  Calculate  a  new  value  of  reference  pressure  drop 
Gf  1  nLx 


(A. 9) 
(A.  10) 

(A.ll) 


AP”ew(ref)  = 


£f=  i«(0A*  k 


(A.  12) 


Step  7:  Adjust  inlet  manifold  pressure  from  new  reference 
pressure  drop 

Pfndj  (ref)  =  AP“w(ref)  +  Pout(ref),  (ref  =  1)  (A.  13) 

P'f  (0  =  Pin 0  -  1)  +  A Pin(i),  0  =  2 to A)  (A.  14) 

Step  8:  Adjust  flow  distribution  factors  from  new  reference 
pressure  drop  and  adjusted  inlet  manifold  pressure 

aadj  ( =  fin^(0  ~  P°utO)  ( ■  =  j  t0  N)  (A.15) 

v  ;  AP"ew(ref)  v  ; 

Step  9:  Calculate  new  set  of  cell  flow  rate  from  adjusted  inlet 
manifold  pressure 

sT  (0  =  [*£j(0  -  ^out(i)]  (i  =  1  to  A)  (A.  16) 
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Step  10:  Check  convergence 


N 

err  =  E 

i=  1 


<?rw(  o  -  sM 

9x{i) 


(A.  17) 


If  ERR  >  e,  update  APr(ref)  and  qx(i)  and  go  to  Step  2.  If 
ERR  <  e,  stop. 
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